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In a previous paper pL a method of comparing the volumes of thermal- 

ized regions in eternally inflating universe was introduced. In this paper, we 
investigate the dependence of the results obtained through that method on 
the choice of the time variable and factor ordering in the diffusion equation 
that describes the evolution of eternally inflating universes. It is shown, both 
analytically and numerically, that the variation of the results due to factor 
ordering ambiguity inherent in the model is of the same order as their varia- 
tion due to the choice of the time variable. Therefore, the results are, within 

their accuracy, free of the spurious dependence on the time parametrization. 
98.80.Hw 

I. INTRODUCTION 

The parameters we call "constants of Nature" can take different values in different parts 
of the universe and in different disconnected universes in the ensemble described by the 
cosmological wave function @||. The probability distribution for the constants can be 
determined with the aid of the "principle of mediocrity" [f|,f|, which asserts that we are 
"typical" among the civilizations inhabiting this ensemble. The resulting probabilities are 
then proportional to the physical volumes occupied by thermalized regions with given values 
of the constants ||, and the preferred values tend to be the ones that give the largest amount 
of inflation. 

This prescription encounters a difficulty when applied to models where inflation is "eter- 
nal". In such models, the universe consists of a number of isolated thermalized regions 
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embedded in a still-inflating background. New thermalized regions are constantly being 
formed, but the inflating domains that separate them expand so fast that the universe never 
thermalizes completely [[F|,[8] . In an eternally inflating universe, the thermalized volumes vi^ 
diverge as t — ► oo. (Here, the index j labels different types of thermalized regions which 
have different values of the constants of Nature.) If one simply introduces a time cutoff by 
including only parts of the volumes that thermalized prior to some moment of time t c , then 
one finds that the ratio 

hm 



is extremely sensitive to the choice of the time coordinate t. For example, cutoffs at a fixed 
proper time and at a fixed scale factor give drastically different results ||. An alternative 
procedure is to introduce a cutoff at the time t^' when all but a small fraction e of the 
co-moving volume destined to thermalize into regions of type j has thermalized. The value 
of e is taken to be the same for all types of thermalized regions, and for all universes in the 
ensemble, but the corresponding cutoff times t%' are generally different. The limit e — > 
is taken after calculating the probability distribution for the constants. It was shown in |T]| 
that the resulting probabilities are rather insensitive to the choice of time parametrization. 

Although this can certainly be regarded as progress, the situation is still not completely 
satisfactory. What do we make of the residual, "mild" dependence of the probabilities on the 
choice of the time coordinate tl This dependence is particularly worrisome when one tries to 
compare probabilities for different universes in the quantum-cosmological ensemble. Using 
the e-procedure, one finds |]] that the probability distribution has an infinitely sharp peak 
at the highest value of the ratio 7/7, where —7 and 7 are the smallest eigenvalues of the 
diffusion equation for the probability distribution of the inflaton field (p, with normalization 
to the physical and coordinate volume, respectively. This sharp prediction does not mix well 
even with a "mild" ambiguity in the value of 7/7. 

An additional source of ambiguity in predictions of eternal inflation is in the form of 
the diffusion equation itself. The equation was introduced in Refs. [|^[||12|-|15|] to describe 



the effect of quantum fluctuations on the evolution of the inflaton field. It gives accurate 
results for sufficiently flat inflaton potentials V (p), provided that the magnitude of V (p) 
is well below the Planck energy density, V ((p) <C 1 . [Here and below we use Planck units 
in which h — c — G — 1.] By its nature, however, the equation is an approximation to 
quantum field theory in a curved spacetime, or, if gravity is to be adequately included, to 
the Wheeler-DeWitt equation for the cosmological wave function. This approximate nature 
of the diffusion equation manifests itself, in particular, in the ambiguity in the ordering of 
factors 8/ dip and D (p>), where D (p) is the field-dependent diffusion coefficient. 

In the present paper we shall analyze the uncertainties in predictions of the principle of 
mediocrity resulting both from the choice of the time parameter and of the factor-ordering 
in the diffusion equation. It will be shown that the uncertainties of the probabilities due 
to these two ambiguities are of the same order of magnitude. Since the factor-ordering 
ambiguity is inherent in the diffusion approximation, the corresponding uncertainty can be 
regarded as the bound on the accuracy allowed by the model. Our result is then that, within 
this accuracy, the probabilities are independent of time parametrization. This is a hopeful 
sign, since it suggests that a complete independence of time parametrization can be achieved 
in a more fundamental approach. 

The paper is organized as follows. In the next section, we introduce the necessary 
formalism. Sections |IT| and |V| present our analytic and numerical results. Conclusions and 
a discussion of some open issues follow in section |V[ 

II. THE DIFFUSION EQUATION 

We shall consider "new" inflation with a potential V (p) of the form illustrated in Fig. [l|. 
The potential has two minima and the values (p* and p^ near the minima correspond to 
the end of inflation. We will be interested in the relative probability 7 ? ( 2 )/7 :, ( 1 ) for the two 
minima. In this section, we shall first introduce the parameters a and j3, representing the 
freedom of time parametrization and factor ordering. We shall then bring the diffusion equa- 



3 



tion to a convenient self-adjoint form, and finally derive a general expression for J>( 2 >/'p( 1 ). 



aK(9) 




FIG. 1. Inflaton potential for a new inflationary scenario. 



A. Parametrization of choices of time and factor ordering 



The evolution of the inflaton field (p during inflation can be described by a probability 
distribution P (ip, t) dip which is interpreted, up to a normalization, as the (either co-moving 
or physical) volume of regions with a particular value of <p in the interval dip. We shall use 
the notation P (ip, t) for the co-moving and P ((p, t) for the physical volume probability. The 
latter satisfies the diffusion equation which can be written as 



is the Hubble constant, V (ip) is the inflaton potential, D (<p) is the diffusion coefficient, 



d t P = d v D (<p) 



l/2+ P d v (D (tpf 2 -? P) - v {#) P] + 3H ( V T P 



(1) 



where 




D(<P) 
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Hto) 



a+2 



(2) 



and 



(3) 



is the average "drift velocity" corresponding to the slow roll. Eq. ([I]) is accurate provided 



that the condition of slow roll, 
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\H'{(p)\ «27riJ(^), (4) 

is satisfied. 

The initial distribution P (ip, t = 0) can be derived from the cosmological wave function 
of the nucleated universe For the "tunneling" wave function, 

P tet = )<xexp(-^). (5) 

Essentially the same results are obtained by taking any bell-shaped initial distribution 
peaked around the maximum of the inflaton potential V (ip). 

The solution is subject to the following boundary conditions at thermalization points 

(1,2) 

(p = (pi ■ 



D(<p) 1/z+p d ip {D(<p) 1/ '- p P) =0. (6) 

The boundary condition (|6|) means that diffusion is constrained to vanish at points (pi ' 
which are found from the "thermalization condition" , 

H'(^)\^2nH(^). (7) 

This equality signifies the breakdown of the slow roll condition (Q). Since (0) is an order-of- 
magnitude relation, the exact values of (pi depend on the choice of the constant of order 
1 in d?]). Although this introduces an ambiguity in our model, we note that diffusion, which 
represents quantum fluctuations, is already negligibly small in the regions dominated by the 
slow roll, and the ambiguity in the choice of points (pi at which diffusion is constrained 



to exactly vanish does not significantly influence the solution of the diffusion equation |T7 . 
The "no-diffusion" boundary condition (||) was introduced in 0. 

In equation (|1|), different choices of time variable and factor ordering are represented by 
choices of the parameters a and (3. The parameter a is equal to the physical dimension of 
the time variable t, which is related to the proper time r by 

dt = H {<p) l - a dr (8) 



(so that a = 1 corresponds to the proper time parametrization, t = r, and a = corresponds 
to using the scale factor as a time variable). The choice of the factor ordering is parametrized 
by (3 which is defined so that (3 = gives the so called Stratonovich factor ordering PJl"2|j . A 
similar parametrization of factor ordering in the diffusion equation ([]]) was used by Salopek 
and Bond . 



The equation for the coordinate- volume distribution P {f,t) is identical to ([!]), except 
for the absence of the expansion term: 

d t P = d v [D {p) 1/2+(3 (p ( V ) 1/2 - p P) - v (<p) P] (9) 

More generally, one could consider a time variable related to the proper time r by 

dt = T (if) dr, (10) 

with T (ip) an arbitrary function; the choice (§) corresponds to T {if) = H (v?) 1_Q . With a 
general time function, the diffusion coefficient would change to 

and the drift velocity would become 

"M=-sfb- (l2) 

Also, the factor ordering in Eqs. (|l]),(P) could be generalized to insertion of an arbitrary 
function h {if) which would change the diffusion term in those equations to 

D {ff^ d v (D {f) 1/2 ^ P)^D {ff 2 h M d v (d {ff 2 . (13) 

The factor ordering of Eqs. ([[]), (|9]) corresponds to the choice h {if) = D {f) . 

We will not consider these possibilities here, since the physical reasons for the "correct" 
choice of functions T {ip) and h {<p) are not clear. In the following discussion, we shall use 
exclusively the parameters a and (3 to explore the ambiguities related to these choices. 
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B. Self-adjoint form of the diffusion equation 



The late-time behavior of the distribution functions described by Eqs. ([!]), (HJ) is deter- 
mined by the eigenfunctions 



P (p, t) = e"^V {1) (<p) , P (<P, t) = e^V (1) (<p) 



(14) 



having highest eigenvalues 7 and —7. The equations for ip (<p) and ip (<p) are stationary 
forms of©,©, 



D (<p) 1/2+p d v (D (tp) 1 '*-" P)-v(<p)P + 3H (<p) a P = ^P 



,1/2-/3 



a- 



D ( V ) l ' 2+f5 d 9 (D (<^) 1/2 - /3 P)-v (<p) P] = -lP 



(15) 
(16) 



The sign convention in the definition of 7 and 7 reflects the fact that the highest value of 
7 and the smallest value of 7 are both positive. (If the highest value of 7 is negative, then 
there is no eternal inflation. See also section [III A] below.) 

Equations fln^fllB]) can be transformed to a manifestly self-adjoint form One in- 
troduces a new independent variable z, 



z = I [2D Mr V2 dip, P M = [2D (v)}'^ f (z) 



1/2 



(17) 



and a further substitution 



f(z,t)=g(z,t)D(zf 2 ex P [— — 5 



= g(z,t)exp (J b(z) dz^j , 



(18) 



where we have defined 



,y \ H' z (5D' Z 

b(z) = -7T— I + - — . 

V ; H 3 2 D 



Equation (H) then leads to the following equation for g (z,t): 

dt9 = \d zz g+{-\v z -\b 2 )g. 



(19) 



(20) 



This form of the equation was derived in [ITJJ for the case of a = (3 = 



7 



Under the same transformations, equation ([!]) for the physical-volume distribution be- 
comes 



d t g=\d zz g+(-\b' z -\b 2 + ?,H a ) tj . 



The transformed versions of Eqs. flI5D,flIED can be written as: 



(21) 



V 

2 y 



l -b' z + ~6 2 - 3H a ) ,) 



-79, 
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(22) 
(23) 



where primes denote differentiation with respect to z. The boundary conditions (ffl) become 



d z 



0. 



(24) 



(1 2) (1 2) 

and they are to be imposed at points zl ' corresponding to the thermalization points (pi ' 
of ©. 

Equations ([Z2"D , (|2"3"D have the form of a stationary Schrodinger equation for a one- 
dimensional motion in a potential U (z), 



with 



-g" + U(z)g = Eg, 



JJ (z) = -b' + -b 2 - 3H a 
K ' 2 2 



for the physical- volume distribution and 



U (z) = -b' + -b 2 
K J 2 2 



for the coordinate-volume distribution. With boundary conditions 
(—1/2) d zz appearing in the left hand side of Q2"5D is non-negative if 



(25) 



(26) 



(27) 



the operator 



/ 9{z) {^2 dzz9 ) dz 



-g {zf ^ (2vr + H 2 f3 (a + 2] 



n4 2) ^ 

+ - [ [g'(z)} 2 dz>0 (28) 
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for any function g (z) . This holds for a potential of Fig. i] since V (<4 2) ) < and 
V (dtA > and so H' z (zf ] ) < and H' z (zi 1] ) > 0. Therefore, the eigenvalues E of 



( Pq) are bounded from below by the minimum values of the potentials ( pq) and and we 
denote —7 and 7 the lowest eigenvalues of the respective operators, as defined by (f22]),(f23l), 
with eigenfunctions ijj^ and t// 1 -*. Exact solutions for ip^ and ijj^ can be found for special 
cases of the inflaton potential V (<p) (and special values of the parameters a, f3) for which 
the Schrodinger equation is exactly solvable (see Appendix A). 

Since the boundary conditions ( |24|) are homogeneous, the operator (— 1/2) d zz + U (z) 
in (|25|) is self-adjoint, and its eigenfunctions are orthogonal with respect to the usual scalar 
product, 



( gi (z),g 2 (z)) = J gi (z)g 2 (z)dz. (29) 
In the original variables, this scalar product becomes 



(P, {fp) , P 2 &)) = JP 1 P 2 ( V ) \~^r e*P dtp. (30) 



2D M 

D{<pf —\ Hto'' 
Using the scalar product (|30|) , the solution of the time- dependent equation (|9]) with a 
given initial distribution P (p, t = 0) = Pq (p) can be expressed through the orthonormal 
eigenfunctions ip^, tp^ 3 \ ... with eigenvalues 7(1) = —7, 7(2), 7(3), as 

00 

P (p, t) = £ C n exp (l(n)t) 0) , (31) 

n=l 

with the coefficients C n given by 

C n = (P to^ {n) ^)), (32) 

and similarly for P (<p,t). Note that co-moving and physical volumes coincide at t — 0, and 
thus P{cp, °) = p ( ( f, 0). 



9 



C. The ratio of physical volumes, r 



The quantities of interest to us are the volumes V* of the thermalization hypersurfaces 
if = (pi . To express these volumes in terms of the eigenfunctions (1^), we first rewrite the 
diffusion equations as continuity equations: 



d t P = -d^J + ?,H a P, 



d t P = -8^ J. 



(33) 
(34) 



Here, the fluxes J, J are related to the distribution functions P, P by 



J=-D ( V ) 1,2+P d v (D (tpf 2 -? P)+v (<p) P. 



(35) 



From the continuity equation (^) it follows that the loss of co-moving volume in inflating 
regions is compensated by the corresponding growth in the volume of the thermalized regions. 
The rate of this growth is given by the flux through the thermalization points: 



0_ 

Of 



(2) 



[ P( V ,t)d V =(-J)\f ) . 



(36) 



„(!) 



(1 2) 

Therefore, the co-moving volume V} C L thermalized in a specific minimum of the inflaton 



(c)* 

potential throughout a given time range is 

*2 



V, 



(1,2) 
(c)* 



J{^ 2 \t) dt 



(37) 



(the absolute value is needed to cancel the negative sign of the flux through the leftmost 
thermalization point). We use a similar formula for the thermalized physical volume: 

*2 



y(i,2) 



According to the method proposed in the infinite thermalization volumes 



(38) 



,(1,2) 



jfaWt) dt 



,(1,2) 



v (<pW) P (<pM,t) dt 



(39) 
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are to be cut off at times t^' 2 ' which are determined from the condition that only a fraction 
e of the corresponding co-moving volume is left unthermalized at those times, 



J(^' 2 \t) dt 



(1,2) 



j(^\t) dt 



(40) 



Here, p± and P2 are fractions of the total co-moving volume (averaged over an ensemble of 
universes) that will eventually thermalize in the first and the second minima of the inflaton 
potential, respectively 

The integrals in (^9f) , (^0|) are dominated by the ground state eigenf unctions ( |14|) of 
(|l5|),(|T6l) which we have denoted by (tp) and ip^ {(f). Using (|3l"l) for P ((p*\t\, we 
obtain 



J {^ 2 \ t) dt= J v (^)) p t) dt « c ^(^' 2) ) v ^(i, 2) ) ^(D ^ 



(41) 



After solving (|p]),(gT|) for 4 1,2) , 

one can calculate from ( |5UD the ratio of the physical 
volumes that thermalize in the two minima of the inflaton potential: 



Vi 



(2) 



r = 



\v 2 \ip2 {\v2\1p2P1 



7/7 



(1) 



\Vl\lpx \\vi\l/jlP2, 



(42) 



where 



v*, a = 0i 1,2) ) , A, a = 0^ 2) ) 



(43) 



(1 2) 

are the values of the ground state eigenfunctions taken at thermalization points y?* , and 



Vl,2 



47T 



J 1 . 2 ) 



(44) 



are "drift velocities" at those points. 



D. Ambiguities in r 



The form of the diffusion equations ([]]), (|5|) depends on the parameters a and /3, and 



therefore the quantities 7, 7, ipj, ifjj, pj appearing in Eq. (|42 ) for the volume ratio r all 
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depend on a and (3. This dependence, which reflects the sensitivity of r to the choice of the 
time parametrization and of the factor ordering, is the main focus of the present paper and 
will be analyzed in detail in the following two sections. 

Another source of uncertainty is the choice of the initial distribution P (ip, 0) = Pq (ip) 
for the calculation of pj. If we choose a Gaussian distribution peaked at the maximum of 
the potential V (ip) and having a width w much smaller than the characteristic width of the 
maximum, then the values of pj should not be very sensitive to w. But a weak dependence 
of pj on w certainly exists and should be addressed as a matter of principle. Some numerical 
results on this w-dependence will be presented in Section |V[ 

The problem of initial distribution can be resolved by invoking quantum cosmology and 
using a distribution, like Eq. ([|), obtained from the cosmological wave function. However, it 
has been argued |19] that probability is an approximate concept in quantum cosmology and 
can be defined only within the semiclassical approximation. The accuracy of this approxi- 
mation for a nucleating universe is characterized by S^ 1 = Hq/tt, where S is the tunneling 
action and H is the value of H at the maximum of V (if). Hence, it may be impossible to 
reduce the uncertainty in P to below O (Hi). 

Yet another uncertainty in Q4"2"D, this time rather benign, is related to the already men- 
tioned choice of the constant of order 1 in (0) which affects the exact values of thermaliza- 
tion points <p^' 2 \ A different choice will change the behavior of the solution tjj^ (ip) near 
ip = y?l 1 ' 2 ' ) (the derivative ■0' (ip) changes sign at y?! 1,2 '* [§]). However, since the diffusion term 
in (|35|) is negligible at those points, one can employ the asymptotic formulae [|J 



^ (1)( ^ 5 f^ eXP 



r* H(ip)d<p _ rv H 1 -* (g) dip 
71 J U H'(tp) + ni J vo W (ip) 



m , s H(ip) 

* { ^ w °1p§) exp 



-47T7 



v H x - a (ip) dip 
H'(ip) 



(45) 



(46) 



which are accurate in the deterministic slow roll region up to the neglected diffusion term. 
Here, <po is some point in the region where diffusion is negligible, 

PjTT 

H 2 fa) « to) ■ (47) 
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With the help of (|45| ) - (|46|) , one can show that a change in thermalization points <£>i — > 
changes the ratio fl42|) by the ratio of the additional volume expansion factors gained before 
thermalization, 



exp(-127r/^ ™^ 



(48) 



To obtain unambiguous relative probabilities, one could compare the volumes of constant- 
temperature hypersurfaces with the same value of T in different types of regions. Then one 
would have to multiply V*** by the additional expansion factors up to the chosen temperature, 

(1 2) 

and the dependence on the precise values of (pi would disappear. 

III. ANALYTIC ESTIMATES 
A. Estimate of eigenvalues 

Perhaps the most important parameters entering (^), in terms of their effect on the 
magnitude of r, are the ground state eigenvalues of Eqs. fllBT) , (p~6|) . These eigenvalues 7, 7 can 
be estimated if the inflaton potential is sufficiently flat and smooth near its maximum. The 
estimate is based on expansion of the effective potential U (z) of the Schrodinger equation 
( p5|) around its minimum up to terms quadratic in z. 

We shall assume for simplicity that the inflaton potential is symmetric around its maxi- 
mum at <p = up to terms of quartic order in (p. Then the expansion of H (<p) around (p = 
has the form 

H M = H + + ^ + ° 0*) ( 49 ) 

with Hq > and H2 < 0. The diffusion coefficient (fj) is also expanded as 

D ((p) =D + + O {(p 4 ) , (50a) 

m +2 (a + 2) m +1 H 2 , . 
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We then use ([H]) to find the derivatives with respect to z at the point z = corresponding 
to <f — 0, 

«9 2 # 



dz 
d 3 H 



dz 3 



2=0 



2=0 



0. 
0. 



dz 2 
d 4 H 



dz 4 



2=0 



2=0 



2D H 2 , 



8D D 2 H 2 + AD 2 H A) 



(51a) 
(51b) 



and substitute them into 
Schrodinger potential, 



6]) and fl27|) to obtain the expansion of the effective 



U(z) = U + ^z 2 + O(z 3 ). 



(52) 



The coefficients Uq and to for the potential (|27|) are given by 



Un 



1 8tt V 2tt / 



rr 

u 2 = H 2H ° 



IT 
2a-2 



167T 2 



2tt 



-5 + 2a + 2(3 (2 + a) + 



(53a) 
(53b) 



and for the potential (^6| 



Uo = U - 3H« = -3H« + \H 2 \ 1 - /3# 2 



2 16tt 2 



AT 


i 


8tt 


-( 






2vr 





r2 a + 2 
2vr 



-5 + 2a + 2/? (2 + a) + 



m 



(54a) 

+ 0(H 4 )). (54b) 



This analysis can be generalized to non-symmetric potentials with H'" (0) 7^ 0; in that case, 
the minimum of U (z) will be shifted from z = 0. 



Assuming that the quadratic expansion fl52|) of the potential U (z) is accurate enough 
up to the classical turning points, we can approximate the eigenvalues of (^5| ) by the corre- 
sponding eigenvalues of the harmonic oscillator, 



E n = U + - + nu, 



(55) 



and the ground state eigenvalues by 



1 = U + -, 



00 



1 = U + -. 



(56) 
(57) 
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Although our boundary conditions are not the same as those for the harmonic oscillator 
eigenvalues (p5[) , they are imposed at points <p* ' deeply within the classically forbidden 
region of the Schrodinger equation and it seems reasonable to assume that a different choice 
of these boundary conditions does not significantly alter the eigenvalues. The validity of the 
estimates fl5E|)-(|57D was confirmed numerically (see Section fV\ ). 

We assume that the potential V (ip) is flat near its maximum and that the maximum 
value V (0) is small in Planck units: 



H 



W"\ 



We also assume that 



2V 



HjHp 



<p=0 



Hi 



in 



V(0) 



<C 1. 



(58) 



(V" 



< l 



(59) 



ip=0 



which is, for instance, the case for the potential 

V{v) = \(p 2 -? 1 2 )' 



(60) 



Under these assumptions, and if the values of a and (3 are not unreasonably large, the 
estimates d5*B|)-(p7|) give: 

\H 2 \ 1 



7 w 3^ 



1 



# 24vr 



\ 



1 + 12a- 



O (Hi) 



(61) 



and 



7 



I -^2 | #0 (-y Hi 



H 4vr 



•>7T 



-5 + 2a + A(3 (2 + a) + + O (h*) 



(62) 



The ratio of the ground state eigenvalues is therefore estimated as 



7 

- w YItx- 

7 -"2 



1 Ho 



1 + 



\ 



1 + 12a- 



m 



Ho 



+ OlH{ 



(63) 



24tt H 

Here, the dependence on a and (3 has been absorbed in O (Hi), except for the a-dependent 
term under the square root. Note that it follows from (|58|),(|6"3"D that 7/7 3> 1. 

To evaluate the a-dependence of the ratio (pBI) resulting from the square root term, we 
have to consider the magnitude of 12Hq/ l-f^l- If 
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H 3 

12^L « i 
\H 2 \ 

then the square root in ( j63| ) can be expanded in powers of Hq/ \H 2 \, which gives 



(64) 



7 in H ® 

— ta I2n- 

7 



Ho 



I Ho I Hn 



12ttH 16tt 



-5 + 6/3 (2 + a) + 



H^Hq 



+ ... 



(65) 



Here we have explicitly written the {Hq) terms, and the ellipsis represents higher-order 
terms. In the opposite limit, 



tt3 

12 wi » x > 

\rLo\ 



(66) 



the a-dependent term in square brackets in (|63|) becomes (aifo \H 2 \ /487r 2 )^ 2 <C H^y/a, and 
we obtain 



1 « 127T H ° 

7 



24vr# n 8tt 



+ ^ -5 + 2a + A/3 (2 + a) + 



H A H 

~hT 



+ ... 



(67) 



We see that, in both cases, the dependence on a and (3 appears only in terms of order O (Hq) 
and higher. 



B. Accuracy and limits of applicability 

The estimates (|6l|)- (|63|) for the eigenvalues of the diffusion equation are valid if the infla- 
ton potential is sufficiently smooth and flat near its maximum, so that the effective potential 
of the Schrodinger equation (p5]) could be approximated as in ( |52j) within a sufficiently large 
region including the classical turning point of (|25|) . This holds if the z A term in the expansion 
of U(z), 

TT ( \ TT , ^ 2 2 , ^4 4 tt $ \H 2 \ 3 Hq^ 1 f 2« H 4 Hq / 2 \\ 

U(z) = U + Y z +-z, U A = i6?r4 — + O (H )), (68) 

is smaller than the quadratic term, (u> 2 /2) z 2 , at the classical turning point z found from 

Uo + y z o =Uo+ 2> Zo = (69) 

which gives the condition 

16 



Gtt 



1 - 



2a H^H 



5m 



« i. 



(70) 
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This inequality holds as long as conditions Q58D and fl59|) are satisfied and a is not very large. 
One can show that the analogous condition for U (z) holds as well. 

Another assumption we made in deriving (|5Bp-(|57|) was uj 2 > 0, uj 2 > 0, meaning that 
the potentials U (z) and U (z) have a minimum (and not a maximum) at z — 0. If (|58| ) -(|59| ) 
hold and a and /3 are of order 1, then uj 2 > also holds, and the condition uj 2 > becomes 

l#2| 



# 3 

12a— ^ > 0, 

\Hn\ 



a > 



12H$- 



(71) 



This condition may be violated for large negative a, but it holds for \a\ ~ 1 if the condition 
(|64]) is satisfied. 

Now we consider the accuracy of the expression fl42|) for the volume ratio r. Since it 
contains the ratio 7/7 in the exponent, and the estimate ( |63D for 7/7 gives an ambiguity 
of order O {Hq), the result of ([EJ) is generally reliable only with logarithmic precision, i.e. 
lnr is determined with a relative accuracy of order O (H^). However, if (Q) is satisfied, the 
a- and /3-dependent terms in ( |65D will be much smaller than 1, and the result of (^) itself 
will be accurate up to O (Hq/ l-f^D- I n this limiting case, we are able to derive a simpler 
approximate formula for the volume ratio r. We notice that under the condition (|64D , the 
region where the quadratic expansion (|52| ) of the potential U (z) is valid, 



-z* > 



a 



4 ~4 



Z < 



2 24 ' 1 1 sJUi 

overlaps with the region where diffusion is negligible 

dH 



(p\ < 



Hp 
\H,V 



(72) 



dip 



\H 2 <p\ , \<p\ > (fo 



Hi 
\Ho\ 



(73) 



Then, the ground state eigenfunction ijj HO (</?) of the harmonic oscillator potential ( p2| ) . 
being a good approximation to the eigenfunction of (|25| ) in the region (72j), can be matched 
with the asymptotic solution in the no-diffusion region, 



^ (if) 



w(v?) 



exp 



-47T7 



1-0 



—dip 



(74) 
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Since the match points (p^ lie within the region ([F^), and the solution ipno (v 9 ) i s symmetric 



in that region, this gives = c^ 2 \ Analogous results are obtained for the physical volume 



distribution, ip^ The formula 
cV\ becomes [ffll 



for the volume ratio r, written in terms of coefficients 



g(2) 



r(l)' 



(2) 



(75) 



With the additional assumption which we verified numerically, 



(76) 



Eq. (|75|) becomes 



K2) 



1 + 



-"o 

l#9.l 



(77) 



where 



Z (l,2) 



exp 



-47T 



±</>0 



(1,2) 



(78) 



are the volume expansion factors during deterministic slow roll. The value of (fo is unim- 



portant as long as it lies within the region (72) where the potential is symmetric. 



IV. NUMERICAL RESULTS 

To check our analytic estimates of the eigenvalues (|6lQ ,(|6^) and to find the dependence 
of the volume ratio (f|2|) on the parameters a and /3, we performed numerical calculations. 
The inflaton potential was chosen as 

V(<p) = V (l-£ + K^ «p(-/^)l . (79) 

where dimensionless parameters k and fi characterize the asymmetry of the potential. We 
considered both the symmetric = 0, and the non-symmetric case (k,{i ^ 0). The 

calculation of the volume ratio ( f£2"D was performed in the non-symmetric case, since it is 
identically equal to 1 for a symmetric potential. 



The conditions d58|)-(|59|) for the potential (|79|) are satisfied if 

V < 1, rf > 1. 



(80) 



A. Technique 



To find the eigenvalues of the stationary equations (p2|) , (|23l) , we used the standard 4th 



5th order Runge-Kutta method |24 |. To facilitate the solution when P (<£>, t) varies greatly 
in order of magnitude, we rewrote the stationary versions of the Eqs. (|33l)-([35"l) in the 
dimensionless variables lnP(<£>, t) and (J / P) as follows: 

d v (In P) = -Svr 2 (p) -(ar + (f + l) # 2 ) a ~ 3 , (81) 
^f^) =7-(^)^(lnP), (82) 

" - N "2 1 ttI Tja-3 



7T+ \2 +1 ) J H ' H ' ( 83 ) 
<9„(T) =3^-7-(i)^(lnP), (84) 

and solved for the smallest values of 7 and —7 to satisfy the boundary conditions, 



J 
P 



(85) 

9* 



where is given by ([H]). The resulting eigenfunctions ift^ 1 ', ip^ were used to calculate 
the values ipip, ^1,2 of (p43|) . 

For the numerical solution of the time-dependent equations (|1|) and (|9|), we have used 
the (slightly modified) unconditionally stable Crank-Nicholson finite difference scheme |2"4j| 
with boundary conditions (f|) and the initial distribution given by a Gaussian, P (if, t = 0) oc 
exp (—{p 2 /w 2 ), with the width u> 77. The solution P (</>, t) was used to obtain the values 

Pl,2 Of ©■ 
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B. Symmetric potentials 

The numerical calculation of the eigenvalues 7, 7 for symmetric potentials ([79]) was 
performed to verify the estimates (|61| ) -(|62| ) . The numerically obtained eigenvalues 7, 7 and 
deviations from the estimates 70, 70 for H = ^J^kVqJ^ = .05 and r\ = 8 as well as the 
ratio 7/7 are summarized in the following tables. The eigenvalues were found with relative 
precision of 10~ 7 . 



a 


-1.0 


-0.5 


0.0 


0.5 


1.0 


7 


0.04982 


0.01114 


0.00249 


5.57e-04 


1.2445e-04 


(7-7o)/7o 


4.94e-04 


9.71e-05 


2.99e-04 


6.9e-04 


8.31e-04 


1 


59.9694 


13.4067 


2.9975 


0.6702 


0.14985 


(7 -To) /To 


7.52e-06 


5.87e-06 


3.05e-07 


2.14e-07 


2.e-07 


7/7 


1203.60 


1203.59 


1203.70 


1203.83 


1203.98 



TABLE I. Eigenvalues for (3 = -1. 



a 


-1.0 


-0.5 


0.0 


0.5 


1.0 


T 


0.04980 


0.01113 


0.002489 


5.565e-04 


1.244e-04 


(t~To)/to 


6.94e-04 


5.94e-04 


4.95e-04 


3.97e-04 


2.97e-04 


T 


59.9694 


13.4068 


2.99751 


0.6702 


0.149851 


(7 -To) /To 


6.69e-06 


2.12e-07 


7.e-16 


l.e-16 


2.e-07 


ill 


1204.08 


1204.07 


1204.17 


1204.31 


1204.45 



TABLE II. Eigenvalues for (3 = 0. 
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OL 


-1.0 


-0.5 


0.0 


0.5 


1.0 


I 


0.04978 


0.01113 


0.002488 


5.5628e-04 


1.243e-04 


(7 -70) ho 


8.93e-04 


1.09e-03 


1.29e-03 


1.49e-03 


1.69e-03 


7 


59.9695 


13.4068 


2.997 


0.6702 


0.149851 


(7-7o)/7o 


5.85e-06 


2.e-07 


3.52e-07 


3.95e-07 


4.17e-07 


7/7 


1204.56 


1204.55 


1204.65 


1204.78 


1204.93 



TABLE III. Eigenvalues for /3 = 1. 



While the eigenvalues themselves vary significantly with a, the ratio 7/7 is very nearly 
constant. One can see that the variance in the eigenvalue ratio 7/7 due to changes in factor 
ordering parameter (3 is comparable to the variance due to changes in the time variable 
parameter a. We have performed numerical calculations using other values of H and 77 in 
the ranges H = 0.5 — 0.001 and r\ = 1 — 20 and obtained similar results. 

C. Asymmetric potentials 

Here we present our numerical results for the potential ([791). 

For H = .05, T] — 5, k = 20, and /i = 50, the eigenvalue ratios and the values of the 
volume ratio ( f42"D are summarized in the following tables. 



a\P 


-1.0 


0.0 


1.0 


-1.0 


469.5152 


469.7020 


469.8889 


0.0 


469.5978 


469.7846 


469.9716 


1.0 


469.6986 


469.8855 


470.0726 



TABLE IV. Eigenvalue ratio 7/7. 
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a\P 


-1.0 


0.0 


1.0 


-1.0 


279.1466 


279.2707 


279.3949 


0.0 


278.1129 


278.2344 


278.3560 


1.0 


277.6672 


277.7888 


277.9106 



TABLE V. Logarithm of the volume ratio (||). 
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As the tables show, the relative change in the eigenvalue ratio is ~ 10~ 3 , which agrees 
with the estimate O (Hq) « 2 ■ 10~ 3 . Also, it is clear that the dependence on a is of the 
same order as the dependence on (3. Calculations were performed for other values of the 
parameters with the same conclusions. 

The form (|79[) of the inflaton potential was chosen to allow for analytic estimates of 
section [III A| , namely, the assumption that the third derivative H'" (0) vanishes holds for 
this potential. We have performed numerical calculations for a potential with nonvanishing 
H'" (0) and obtained similar results. 

To verify the assumption (|76|) used in our derivation of (j77|), we looked at the dependence 



of In {px/p-i) on a and (3. Our results suggest that the relative variance of In (px/p?) with a 
and /3 is also of the order Hq. However, in our case of a "flat top" potential, the symmetry 
of the potential in the diffusive region leads to pi ~ P2, so that 

^ln^«l, (86) 

7 P2 

and the volume ratio (|77|) is virtually unaffected by the dependence of In (pi/pz) on a and 
(3. For potentials with asymmetry in the diffusive region, the value of In ^ was of the order 
1, however its dependence on a and (3 remained small (O(Hq)). 

Another possible source of uncertainty discussed above was the choice of the initial 
distribution. We performed calculations of the volume ratio (|42| ) for a Gaussian initial 
distribution, 

2" 



P ((p) oc exp 



WT] 



with varying width parameters w = 0.0001 — 0.01, and the results varied insignificantly (the 
relative change in lnr was of the order 10~ 5 ). 
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V. CONCLUSIONS 



In this paper we analyzed the ambiguities in assigning probabilities to different types of 
thermalized volumes in an eternally inflating universe. One of the key factors determining 
the relative probabilities is the volume ratio r, given by Eq. (0). Our results are most 
easily formulated in terms of the quantity y = lnr. 

We introduced a parameter a representing the ambiguity due to the choice of the time 
variable t, and parameter f3 representing the factor ordering ambiguity in the diffusion 
equation ([!]). Our main result, obtained both analytically and numerically, is that variation 
of either of these parameters introduces a variation 

Sy/y ~ Vo, (87) 

where Vq ~ Hq is the highest value of the inflaton potential V (if), and H is the expansion 
rate of the universe at V = Vq. 

Since the factor ordering ambiguity is inherent in the diffusion approximation, Eq. ( |87|) 
gives a bound on the accuracy of the predictions of the model. Moreover, since variations 
of y arising from varying a and (3 are of the same order, we conclude that within that 
accuracy, the results are independent of the choice of time variable. This is an intriguing 
result, since it suggests that in a more fundamental approach, based e.g. on the Wheeler-De 
Witt equation, probabilities could be manifestly independent of time parametrization. 

Another source of uncertainty is the choice of the initial distribution for the calculation 
of pj in Eq. fl42|). Our results indicate that the corresponding variation of y is even smaller 
than fl87|). One could try to avoid this uncertainty by using an initial distribution, like Eq. 
(H), derived from the cosmological wave function. However, probability is an approximate 
concept in quantum cosmology, and unitarity holds only with the accuracy of the semiclas- 



sical approximation For a nucleating universe, this accuracy is of order Vq (see Sec. 

HDD . 



Throughout the paper we have assumed inflation of the "new" type with a potential 
well below the Planck scale, V (if) <C 1. In this case, Sy/y <C 1. As the maximum value of 
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the potential Vq is increased, the accuracy gradually deteriorates, and errors become 0(1) 
when it reaches the Planck scale, Vo ~ 1. This is expected to happen in models of "chaotic" 
inflation, where the probability distribution P (ip, t) is concentrated near the Planck values 



of the potential [20,21 



The relative uncertainty in the volume ratio r itself is typically greater than Q87|), due to 

the presence of the large exponent, 7/7 ~ Vq] \V "\, in Eq. ((42j): 

frr V 2 

T = s »~m- (88) 

7 I v o I 

This is small if V 2 / \ V "\ <C 1, in which case r is accurately given by the simple formula 
(|77D . Otherwise, the uncertainty in r is large. We note, however, that the probabilities of 
thermalization into different minima of V (<£>) are expected to be vastly different, so that 
the corresponding values of y are large, \y\ 3> 1, and the resulting volume ratios are either 
very large (r >>> 1) or very small (r <<< 1). These strong inequalities are unaffected by 
the uncertainty (0). It can affect only rare borderline cases when the two probabilities are 
nearly equal. It appears that we have to accept this as a genuine uncertainty of the problem 
and make predictions only in cases where one minimum is much more probable than the 
other. 

In this paper, we considered exclusively the problem of finding the relative probabilities 
for thermalization into different minima of the inflaton potential in a single universe |30| . 



The same method can also be applied || to an ensemble of disconnected, eternally inflating 
universes with different potentials V ((f), parametrized by some variable "constants of Na- 
ture" {Xj}. (The set of allowed values of {Xj} may be either discrete or continuous.) The 
volume of thermalized regions in a given universe would then depend on the cutoff parameter 



e as 



V, oc e~^, (89) 

where the eigenvalues 7, 7 pertain to the diffusion equation with the potential V ((p) in that 
universe. In the limit e — >• 0, only universes with maximum value of 7/7 will have a non-zero 
probability |22 . 
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It is possible that the condition 

7/7 = max (90) 

selects a unique set of {Xj}- Then the potential V (<£>) is fixed, and the only remaining 
problem is to find the relative probabilities for thermalization into different minima of this 
potential. On the other hand, it is conceivable that the maximum of 7/7 is strongly degen- 
erate p3| , so that Eq. (|9~0|) selects a large subset of all allowed values of {Xj}- It is worth 



noting that the class of potentials having the same "ground state" eigenvalue 7 is very wide: 
it can be parametrized by an arbitrary function (see Appendix B). The relative probabilities 
for {\j} within the degenerate subset can be calculated following the same procedure as in 
Ref. |l| and in the present paper. 

NOTES ADDED 



1. After this paper was submitted, we learned about an interesting paper p8[ , in which 
arbitrary time parametrizations fllCf ) for the diffusion equation ([]]) were considered. The 
authors described a transformation of potential V (if) and a corresponding change of time 
variables that give identical physical predictions at late times. In the framework of the 
present paper, that equivalence transformation leaves the Schrodinger equation ( ^5|) and 
boundary conditions (|24l) unchanged, giving the same eigenvalues and eigenfunctions and, 
therefore, identical predictions. The family (|A7|) of exactly solvable potentials was also found 
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2. In a recent preprint p9| , Linde and Mezhlumian suggested a family of alternative 
regularization procedures parameterized by a dimensionless number q. All these procedures 
have the same property of time reparametrization invariance as the one we discussed here, 
indicating that the invariance requirement alone is not sufficient to select a unique regular- 
ization procedure. In this note we shall briefly discuss some additional requirements which 
may fix the parameter q. 
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The alternative regularizations of |29| are identical to ours, except the co-moving volume 
distribution P ((p, t) is replaced by a weighted distribution P q {ip, t) which satisfies a modified 
version of Eq. (|33| ) , 

d t P q + d 9 J g = 3 gi! ^—P q . (91) 

(The value q = corresponds to the unmodified regularization procedure of JTJ.) For q > 0, 
a greater weight is assigned to co-moving regions which have expanded by a greater factor; 
for q < 0, to regions which have expanded by a smaller factor. The latter situation appears 
somewhat unnatural, and we can require that 

q > 0. (92) 

The parameter q should be chosen so that the smallest eigenvalue 7 g of (^TJ) satisfies 

lq > 0, (93) 



ensuring that the integrals in Eq. ( ^0|) are convergent. Using the scale factor time 
parametrization (a = 0) we have 7 g = 7 — 3g, where 7 = 7^=0 takes values in the range 
< 7 < 3, depending on the inflaton potential V (p). If we require that the regularization 
scheme should apply to all possible inflaton potentials, then the condition ( |9"3"D restricts q 
to the range q < 0. Combining this with the condition (|9~2"D , we are left with a single value, 
q = 0. 

Let us finally consider the situation when the potential V (p) in Fig. 1 is symmetric 

in the range of p where diffusion is non-negligible, so that the difference between V* and 

(2) 

V* is due entirely to the regions of deterministic slow roll. We can require that in this 
case the volume ratio should be the same as in the case of non-stochastic, finite inflation, 
r = {z^ I zi 1 ^ , where Z* 1 ^ and Z^ are the corresponding slow-roll expansion factors. 
Once again, this selects q = 0. 

Although none of the conditions we have suggested appears to be mandatory, the above 
discussion does suggest that the regularization scheme with q = has some unique features 
and may therefore be preferred. 
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APPENDIX A: EXACT SOLUTIONS OF THE DIFFUSION EQUATION 



For specific inflaton potentials V ((f), the Schrodinger equation ([250 is exactly solvable. 
We will present several examples of such potentials. 

The widest range of solvable potentials occurs with the choice a = (3 = 0. In this case, 



equations (22) and ( p3"D differ only by a constant shift of the eigenvalue, and it is sufficient 
to consider one equation (^3). The potential (E7|) is given by 



U (z) = l -b> (z) + h? (z) 



(Al) 



where 



b(z) 



-it 



H'(z) 
if 3 (z) 



(A2) 



The simplest case of a solvable Schrodinger equation is that with a constant (or piecewise- 
constant) potential, U (z) = Uq. This is achieved, for instance, if 

-1/2 



b(z) = ±J2U n , Hz) 



7T 



z = z ± 2tcJ2Uq{p 



(A3) 



which corresponds to 



V (ip) = const, Uq = 0. 



(A4a) 
(A4b) 



Taking V (ip) to be of the form ( |A4a| ) and ( |A4b| ) with different Uq on different regions 
of values of ip, one obtains a piecewise-constant potential U (z) = Uq for the Schrodinger 
equation (|23|) . 
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Another exactly solvable case is that of the harmonic oscillator potential, 

U(z) = U + ^z 2 , (A5) 

with 

/ 2 \ -V2 

[/ = ±|, b(z)=±uz, H(z)=(±^- + c) , (A6) 
The upper choice of sign corresponds to the inflaton potential 



vM = v -p(-^ (A7) 

(1 + cexp ( — A^ixujLp)) 



Particular cases of these potentials include the exponential potential ||25|| , V (<p) = 
Vq exp (cup) , and the potentials 

Vfo0 = — , v b) = -zS-r- ( A8 ) 

cosh cup cos z cup 



APPENDIX B: POTENTIALS WITH A GIVEN GROUND STATE EIGENVALUE 

It is always possible to choose a potential for a Schrodinger equation that would possess a 
given ground state eigenvalue, and the freedom of this choice is parametrized by an arbitrary 
function. One can show this using methods of supersymmetric quantum mechanics [J26[ |. We 
start with an arbitrary function / (z) and define the potential U (z; f, Eq) for the Schrodinger 
equation (|25|) by 

U (z; f, E ) = l - exp (-/ (z)) ^ exp / (z) + E = \ {f + f) + E . (Bl) 

It is easily verified that the function g (z) = exp / (z) satisfies (|25|) with eigenvalue Eq. 
Since the function g (z) is everywhere nonzero, it would be the ground state eigenfunction if 
appropriate boundary conditions were satisfied. Generic homogeneous boundary conditions, 

(g' + B(z)g)\ Zj = 0, (B2) 

imposed at some boundary points Zx,2, translate into boundary conditions imposed on / (z): 
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f{z j ) = -B{z j ). (B3) 

Since the boundary conditions restrict the behavior of / (z) only near the boundary points, 
the freedom of choosing a function / (z) is essentially unaffected by the boundary condition 
requirement. We see that for all functions / (z) the corresponding potential U (z\ f, E ) has 
the required ground state eigenvalue Eq. 

It is also always possible to find an inflaton potential V (if) for Eq. (|]) that will lead 
to U (z; f, E ) in the corresponding Schrodinger equation. To do that, one needs to find a 
solution g (z) of the equation (|25| ) with potential ( (BTj) with eigenvalue (where the function 
go (z) does not have to satisfy the boundary conditions). Then one can define b (z) = g' /g 
and express the potential U (z; f, E ) as in (P7|), 

U(z;f,E ) = ^(b' + b 2 ). (B4) 

Using the terminology of SUSY quantum mechanics, the function —b (z) is the "superpoten- 
tial" . By transforming back to the physical variable if, one can obtain the potential V (if) 
corresponding to this Schrodinger equation. 
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